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LEXINGTON 


MASSACHUSETTS 


ABSTRACT 


In  this  report,  we  extend  a  simple  graphical  model  for 
analyzing  radar  homing  interceptor  engagements  [1]  to  include  the 
radome  refraction  error  effect  and  a  refined  miss  distance 
estimator.  The  residual  radome  error  is  treated  as  an  additive 
exponentially  correlated  noise  with  correlation  time  equal  to  the 
interceptor  response  time.  The  miss  distance  is  estimated  as  the 
error  in  predicting  the  target  cross-range  position  at  intercept 
at  the  range  where  the  interceptor  divert  capability  becomes  just 
insufficient  to  correct  the  distance  difference  between  the 
previous  and  the  current  predicted  intercept  points.  The 
analytical  and  numerical  results  show  that  the  prediction  error 
is  a  strong  function  of  the  correlation  time  and  the  effect  of 
correlation  generally  increases  the  achievable  miss  distance 
especially  when  the  interceptor  response  time  is  large. 
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I.  INTRODUCTION 

There  is  considerable  interest  in  evaluating  the  performance 
of  homing  interceptors  using  simple  analysis  techniques  [1,  2], 
The  analytical  results  are  generally  good  enough  for  system-level 
studies  and  could  be  used  as  guides  for  complicated  computer 
simulations. 

In  [1],  a  simple  graphical  model  for  analyzing  radar  homing 
interceptor  engagements  was  presented.  The  model  compares  the 
accuracy  of  the  intercept  point  prediction  and  the  interceptor 
divert  capability  to  estimate  the  final  miss  distance.  At  long 
ranges,  the  interceptor  can  usually  taKe  out  the  prediction 
error  because  of  the  long  time  remaining  until  intercept. 

However,  at  sufficient  close  ranges,  the  interceptor  is  no  longer 
capable  of  fully  correcting  the  prediction  error  and  this 
residual  error  is  considered  as  a  rough  estimate  of  the  miss 
distance.  In  this  model,  the  sensor  error  is  treated  as  a  white 
noise  and  it  is  assumed  that  the  interceptor  responds  to  the 
sensor  prediction  perfectly.  (This  has  been  called  "divine 
guidance. " ) 

In  this  report,  we  extend  the  previous  model  [1]  to  include 
the  boresight  slope  error  effect  and  a  refined  engagement  model. 
The  boresight  slope  error  which  is  believed  to  be  important  to 
the  final  miss  distance  is  modeled  as  a  correlated  noise  additive 
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to  the  measurement  [3,  4].  The  analytical  and  numerical  results 
show  that  the  prediction  error  is  a  strong  function  of  the 
correlation  time  of  the  noise.  Besides,  the  effects  of 
correlation  generally  increase  the  achievable  miss  distance, 
especially  when  the  interceptor  response  time  is  large. 

The  refined  engagement  model  assumes  that  the  sensor  revises 
the  prediction  of  target  position  at  intercept  at  each  pulse  and 
the  interceptor  changes  maneuver  accordingly  to  get  to  the 
predicted  position.  The  miss  distance  is  then  regarded  as  the 
prediction  error  at  the  range  where  the  interceptor  divert 
capability  becomes  insufficient  to  correct  the  distance 
difference  between  the  previous  and  the  current  predicted  target 
positions . 

In  the  next  section,  we  first  summarize  the  basic  model 
described  in  [1]  in  order  to  provide  a  foundation  for  the 
following  discussions.  In  Section  3,  we  describe  the  refined 
sensor  error  model  and  the  resulting  prediction  errors. 

Section  4  describes  the  final  miss  distance  model  and  presents 
some  numerical  results.  Finally,  a  conclusion  is  given.  There 
are  two  appendices  containing  detailed  derivations  of  the 
statistical  formula  used  in  Sections  7  and  4. 


II . 


BASIC  MODEL 


In  this  section,  we  summarize  the  basic  sensor;  interceptor 
and  engagement  models  are  presented  in  [1].  The  homing  sensor 
tracks  a  target  from  acquisition  range  ra  to  a  given  range  r 
and  predicts  ahead  to  intercept  (range=0).  Based  on  a  sequence 
of  measurements  of  the  target  cross-range  position,  z^ 

(i=1,...,  N),  at  range  r^  with  accuracy  and  assuming  a 
linear  target  trajectory  model 

zi  =  0O  +  e1  ri  +  wi  (1) 

where  w^  i=1,...,  N  is  a  sequence  of  uncorrelated  zero-mean 

noise,  the  general  formula  for  the  predicted  position  {  Bg)  and 

A  /  ~7\ 

the  prediction  error  a(  8  )  =  V  Var ( Bg )  have  been  derived.  Among 
various  radar  error  sources,  we  are  particularly  concerned  with 
the  instrumentation  error.  For  this  error,  0i=oe  r^  where 
oq  is  the  radar  angular  accuracy  and  Bg  and  Var  (Bg)  are 
given  by 


(2) 


3 


and 


2 


Lospectively. 

The  interceptor  response  to  a  commanded  acceleration  is 
modeled  as  a  constant  acceleration  preceded  by  a  pure  time 
delay,  xm .  Thus,  the  divert  capability  of  the  interceptor  at 
range  r  is  equal  to 

D  =  1/2  a(r/vc  -  Tm)  (4) 

where  vc  is  the  closing  velocity  and  a  is  the  maneuver 
acceleration  limit. 

During  engagement,  at  an  instantaneous  range  after 
acquisition  at  which  tracking  stops  and  prediction  starts,  the 
interceptor  is  commanded  to  erase  the  prediction  error.  In  the 
early  stage  of  engagement,  the  divert  capability  of  the 
interceptor  is  usually  more  than  enough  to  do  so.  However  as  r 
decreases,  the  divert  capability  decreases  faster  than  the 
prediction  error  and  so  at  some  point  the  divert  is  just  adequate 
to  take  out  the  prediction  error.  Below  this  point,  the  divert 
becomes  insufficient.  The  corresponding  prediction  error  at  this 
point  is  then  considered  as  an  estimate  of  the  final  miss 
distance. 
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III.  A  BETTER  SENSOR  NOISE  MODEL 


In  this  section,  we  calculate  the  predicted  target  position 
at  intercept  and  the  associated  prediction  error  using  a  refined 
sensor  noise  model  which  can  take  the  radome  refraction  error 
into  account.  We  will  follow  the  same  notation  as  used  in  the 
previous  section. 

It  has  been  pointed  out  that  the  radome  refraction  error 
could  cause  instability  of  the  interceptor  [3]  and  should  be 
compensated  by  an  on-board  computer  using  pre-calibrated  look-up 
tahles.  Even  after  compensation,  the  residual  error  is  still 
expected  to  influence  the  performance  of  the  interceptor.  This 
residual  radome  error  is  modeled  as  an  exponent ial ly  correlatel 
noise  in  the  look  angle  domain  [4].  By  proper  choice  of  error 
magnitude  and  correlation  distance,  it  is  possible  to  match  the 
appropriate  boresite  slope  error  distribution.  In  simulation, 
this  correlated  noise  can  be  generated  by  passing  a  white 
sequence  through  a  first  order  Markov  process  filter. 

We  assume  that  noises  additive  to  the  target  cross-range 
position  measurements  due  to  radome  refraction  are  also 
exponentially  correlated  but  in  the  range  or  time  domain.  That 
is 

E  (w.  w.]  =  o.  o,  pP  1  (5a) 

l  g  l  g 


i 


(5b) 


p  =  exp  (-At/Tc)  =  exp  (-1/(prf  •t(;)  ) 

where  At,  equal  to  1/prf,  is  the  pulse  interval  and  t  is  the 

c 

correlation  time.  If  the  interceptor  does  not  maneuver  by 

changing  its  angle  of  attack,  the  Tc  +  °°*  If  the  interceptor  does 

maneuver,  then  a  reasonable  assumption  for  is  that  it  is 

equal  to  the  missile  response  time,  x  .  This  is  because  the 

look  angles  are  likely  to  be  correlated  within  t  due  to  lack  of 

m 

maneuver  change  but  tend  to  vary  significantly  over  a  period 

longer  than  r  due  to  the  new  maneuver. 

3  m 

According  to  this  error  model,  it  is  shown  in  Appendix  I 

A  A 

that  Bq  and  Var(0Q)  are  given  by 

'b'q  =  a  U(£r  £)  f ( £,  j)  -  f(£»  z)  f(£,  j)  }  (6) 


A 

Var ( Bq ) 


»2  (  1-P2)  a  f  (£,  £ ! 


(7) 


where 


f  ( U  ,  V  )  = 


i'  i> 

f  (£» 

£)  -  f2 

<t.  j))'' 

(  N-l 

(U_i_ 

P  Ui+  l  \ 

1  < 

I  H‘ 

1 

■O 

< 

H 

(  i=i 

\  r. 

1 

ri+l  > 

\  r.  r  .  ,  /  t 

i  i  +  I  / 

•  •  •  9 

1  )  ,  £ 

-  (  r  ^  /  •  • 

.,  rN),  z  =  (z, 

.  a-..2)  ^ 
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when  o . 

1 


=  o0  r^.  Here  o0  might  represent  the  magnitude  of  the 
residual  radome  refraction  error  or  the  radar  range-independent 
(instrumentation)  error.  To  be  more  explicit,  vart^)  in  (7)  can 
be  rewritten 


Note  that  when  p=0  (or  ^=0),  the  case  of  uncorrelated  noise, 
the  above  two  equations  reduce  to  (2)  and  (3)  as  expected. 

Figure  1  shows  the  prediction  error  as  a  function  of  the 
instantaneous  range  r  for  several  values  of  using  a  set  of 
nominal  values  for  rfl(5  km),  vc(5  km/s),  prf(IOO)  and  o0(2  mr ) . 
It  is  clear  that  at  long  ranges  the  correlation  degrades  the 
prediction  accuracy  but  starts  to  improve  it  near  intercept. 

This  behavior  can  be  explained  as  follows:  The  effects  of 
correlation  are  twofold.  On  one  hand,  higher  correlation  means 
fewer  independent  samples  and  so  larger  o(^g).  On  the  other 
hand,  with  short  time  to  intercept,  the  higher  correlation 
imposes  strong  proportionality  constraints  among  samples  [5J  and 
reduces  the  dependence  of  ^  on  the  noise.  In  the  limiting  case 
where  p=1  (or  ^c=") ,  for  any  given  w^,  the  subsequent  noise 
samples  are  restricted  to  be 
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I 


i =2 , .  4  ,N 


Clearly, the  noise  does  not  affect  8^  and  therefore  3^  is  always 
equal  to  8Q  and  o(8g)=0.  At  long  ranges,  the  number  of  available 
samples  is  small  and  so  the  first  effect  might  dominate  for 
moderate  values  of  p.  Near  intercept,  many  samples  become 
available.  The  first  effect  is  less  important  and  is  overwhelmed 
by  the  second  effect.  In  fact,  it  can  be  shown  that  for 
r<<vc/prf, 


0 . 
l 


w .  = 
i 


w. 


and  (1)  can  be  written  as 


w . 


z .  =  3  +  (  8  + 

l  o  1  r. 


A 
o  (  3 


(8) 


However,  this  corresponds  to  the  region  after  the  last 
measurement . 

Figures  2  and  3,  in  comparison  to  Figure  1,  show  the  effect 
of  prf  on  the  prediction  error  while  o0  is  kept  constant. 
Generally  speaking,  higher  prf  tends  to  reduce  the  prediction 
error  for  the  case  tc  =  0  but  not  for  the  other  cases. 
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IV.  MISS  DISTANCE  MODEL 


The  intercept  point  predicted  by  the  homing  sensor  changes 
from  pulse  to  the  pulse.  At  the  Nth  pulse,  the  interceptor 
responds  to  the  revised  prediction  by  maneuvering  to  take  out  the 
difference  between  ®g(N)  and  ^(N-l)  within  the  remaining  time 
until  intercept.  In  other  words,  the  interceptor  is  commanded  to 
chase  the  predicted  intercept  point.  During  the  engagement 
process,  the  instantaneous  miss  distance  becomes  smaller  and 
smaller  as  the  interceptor  approaches  the  target.  At  some  point 
the  divert  capability  of  the  interceptor  becomes  insufficient  to 
respond  to  the  chanqe  in  predicted  target  position.  At  this 
critical  point,  the  prediction  error  can  no  longer  be  reduced  and 
it  is  considered  as  an  estimate  of  the  final  miss  distance. 

The  one-pulse  position  change  at  Nth  pulse  is  a  random 

A  A 

variable.  Let  AN= 3^{ N ) - ( N-l ) .  Then,  the  variance  of  as 
shown  in  Appendix  2,  is  given  by 


Var(AN)  =  var(§0(N-1 ) ) 


A 

-  var ( U0(N) 


(9) 


and  it  can  be  computed  using  the  formula  for  Var ( 3^ )  in  Eq.  (3). 

Statistically,  the  one-o  value  of  A^  is  considered  as  the 

distance  which  should  be  taken  out  by  the  interceptor  after  the 

Nth  pulse  or  range  r=r  -(N-1)  v  /prf. 

a  c 
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An  illustration  of  the  engagement  model  is  shown  in 


Fig.  4  where  the  prediction  error  of  the  sensor  [o(3g),  Eq.  (3)] , 
the  divert  capability  of  the  interceptor  [D,  Eq.  (4)1  and  the 
one-pulse  position  change  [o(^),  Eq.  (9)]  are  plotted  as 
functions  of  the  ranqe  r.  At  very  short  range,  °(AN)  is  larger 
than  o(8q).  This  is  because  the  relative  change  in  between 
pulses  is  significant.  As  long  as  the  time  between  pulses  is 
much  shorter  than  the  missile  response  time,  this  effect  does  not 
influence  the  miss  distance.  The  critical  point  where  the  divert 
and  position  change  curves  cross  over  each  other  is  designated  as 
"A".  The  corresponding  miss  distance,  M,  defined  as  the  predic¬ 
tion  error  at  this  time  is  .4  meter. 

Based  on  the  above  described  miss  distance  model,  sample 
numerical  results  are  presented  in  the  following.  The  nominal 
operating  point  of  the  engagement  is  assumed  to  be  r  =  5  km, 

ci 

vc=5  km/s,  <Jg=  2  mr ,  prf  =  100,  t^=  40  ms  and  a  =  40  q.  The 
miss  distance  variations  corresponding  to  certain  excursions  of 
og,  t  ,  a,  and  prf  are  shown  in  Figures  5-8,  respectively.  Two 
cases  are  considered  in  each  figure,  uncorrelated  noise  and 
correlated  noise.  In  the  latter  case,  the  correlation  time  is 
set  equal  to  the  missile  response  time,  i.e. 

Two  observations  are  very  distinct  from  these  figures. 
Firstly,  the  miss  distance  increases  when  ag  becomes  larger,  when 
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t  becomes  lonqer,  when  a  becomes  smaller  or  when  prf  becomes 
m 

lower.  The  chanqes  are  more  prominent  for  the  case  of  correlated 

noise.  Secondly,  the  effect  of  correlation  deqrades  the 

achievable  miss  distance  in  all  cases  under  consideration.  The 

deqrees  of  deqradation  vary  but  they  are  particularly  siqnificant 

if  t  is  larqe  or  if  a  is  small, 
m 
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Fig.  5.  Miss  distances  for  various  magnitudes  of  the  radome 
residual  error. 
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Fig.  6.  Miss  distance  for  various  interceptor  response  times. 
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Fig.  7.  Miss  distances  for  various  interceptor  acceleration 
limits. 


V. 


CONCLUSION 


We  have  modified  the  basic  model  of  a  radar  hominq 
enqaqement  presented  in  [1]  to  include  the  boresiqht  error  effect 
and  a  refined  miss  distance  model.  The  boresiqht  error  is 
treated  as  an  exponentially  correlated  random  noise  additive  to 
the  tarqet  position  measurement.  The  final  miss  distance  is 
estimated  as  the  prediction  error  at  the  ranqe  where  the  divert 
capability  of  the  interceptor  is  just  sufficient  to  erase  the 
chanqe  in  predicted  intercept  points  from  the  previous  to  the 
current  predictions.  Numerical  results  show  that  the  effect  of 
correlation  qenerally  increases  the  achievable  miss  distance, 
particularly  when  the  interceptor  response  time  is  larqe  or  when 
the  acceleration  limit  is  small.  Finally,  we  would  caution  the 
reader  that  even  with  some  modifications  on  the  basic  model,  the 
model  described  in  this  report  is  still  very  simplified  and 
should  be  used  to  qive  a  lower  bound  on  achievable  miss  distance. 
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APPENDIX  A:  Derivation  of  Prediction  Error  Formula 

■ 

i  In  this  appendix,  we  will  derive  Eqs.  (6)  and  (7)  of 

[  Section  3.  They  are  actually  obtained  as  a  special  case  of  a 

f  rather  general  approach.  Given  a  sequence  of  measurements  z^ 


Let  W  denote  the  covariance  matrix  of  the  noise  sequence. 


i  .e . , 

W  =  E{w  wS 

Then,  according  to  the  well-known  Gauss-Markov  theorem  [6],  the 

A 

best  linear  unbiased  estimate  of  ji,  denoted  by  J,  is  given  as 


i  =  (xfc  w-1  X)"1  xfc  w"1  z 

A  . 

and  the  covariance  of  J5  is 
cov(B)  =  (Xfc  w'1  X)'1 


(A4) 


(A5) 


(A4)  and  (A5)  can  be  simplified  if  W  possesses  the  special 
structure  that 


w 


N-l 


°1  °2  P 


o1  o2  p  .  .  .  •  °N  P 


2  ,  N-2 

. °2  °N  ° 


0.  0. 


N-l 


N-2 

P 


I 


(A6) 


as  described  in  Section  3.  It  is  easy  to  see  that  W  in  this  case 
can  be  decomposed  into 
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w 


DVD 


(A7) 


where 

D  =  diag  .  .  ,  °N ) 

and 

p”-1  \ 

N-2 

P 

1 

By  direct  substitutions  and  using  the  fact  that  [7] 


/ 


l-p 


-P 


o 


2 


-P  1+p  -p 

2 

i+r 

O 


i+p 

-p  i 


both  £  and  cov(£)  can  be  expressed  explicitly  without 

A  A 

inversions.  Specifically,  and  var(3g)  are  given  by 


(A8) 


(A9) 


(A10) 


the  matrix 
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Sq  =  a  {f(r,  r)  f(z,  j)  -  f(r,  z)  f(r,  i) \  (All) 

var ( )  =  (1-p2)  «  f(£,  r)  (A12) 


where 


a  =  {f(i»  2>  f(£/  £)  "  f2<£»  jj  1  1 


f(U,  V) 


(1-p2) 


V 


N 

2 


0 


N 


For  the  case  that  =  oQ  with  o0  being  a  constant  (All)  and 
(A12),  reduce  to  (6)  and  (7)  of  Section  3. 
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APPENDIX  B:  Derivation  of  Formula  for  Intercept  Point 
Prediction  Change  from  Pulse-to-Pulse 


In  this  appendix,  we  will  derive  Eq.  (9)  of  Section  4  using 
the  observation  model  described  by  (A2)  and  (A6).  Two  subscripts 

"N"  and  "N-1"  will  be  used  to  denote  quantities  constructed  from 

A  A 

N  and  N-1  measurements,  respectively.  For  example,  and  J^j_i 
denote  the  estimates  of  at  the  Nth  pulse  and  the  (N-I)th 

pulse.  During  the  derivation,  three  important  relationships  are 

A 

obtained  in  order:  (1)  cov  )  can  be  computed  recursively,  (2) 


A 


A 


cov  {j*N,  )  =  cov  (j»N)  and  (3)  cov  (B_N  mn_1,  -  vmn_x 


/N  /N 

B..  ,  )  =  cov  ( B..  .  )  - 


A 

c°v  (_Bn  )  . 


From  ( A5 ) ,  we  have 

cov  (  0., )  =  ( X  *  W  “1  XM)_1 


-N' 


N 


N  ' 


(31  ) 


Both  X„  and  W,, 
N  N 


-1 


can  be  partitioned  as  follows: 


X  = 
N 


( B2 ) 
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0 


where 


t 

P 


a.  v 


/  o  O  \ 


S  = 


T 


S  /  0 


2 

N-l 


(B4a ) 


(B4b) 


(B4c) 


(0, . . .0, 


(1-°  >  °N  Vl 


(B4d) 


g  =  [ (1-p2)  O2  1  1  (B4e) 

Substituting  { B 2 )  and  (B3)  into  (B1)  and  carrying  out  the  matrix 
manipulation  yields 


26 


■  <XN-.  WN-1  XN-t  +  XS-1  T  XN- 1+  E  Vl* 

XN-1  1  EC  +  E  9  Et,_1  <B5> 

By  direct  matrix  multiplications  and  summation,  it  can  be  shown 
that  the  summation  of  the  last  four  terms  within  parentheses  of 
(85)  is  equal  to 

Tv  %  ~n 


where 


/ 


% 


N-l 


P  r 


N-l  N 


\  °N-1 


V 


(B6) 


Therefore 


A 

c°v( SN) 


(XN-1  WN-1  XN-1  +  “  2  3N  2m* 

1-p 


(XN-1  WN-1  XN-1)  "  °N 


=  cov(iN-i}  "  QN 


(B7a) 
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I 


where 


(1-p2)-1  cov^)  gNt  COV^) 
1  +  (1-P2)"1  qNt  COV(^J_1)  ^ 


(B7b) 


In  the  above,  the  well-known  matrix  inversion  lemma  [7]  has  been 
used . 

A  A 

Next,  we  will  show  that  the  covariance  between  and  0N_1 , 
cov(_3n,  £n_1),  is  equal  to  cov(_$N).  Let 


%  -  <XN  Wn'  V'1  XN  V 


( B8 ) 


Then  the  following  relationships  are  immediate: 

% -  cov(^>  4 

-N  =  UN  -N 


UN  XN  "  1  • 


Using  these  equations,  we  have 
A  A  A  r  A 


cov(0N,  1N-1)  =  e{(Jn_1  -  £)  (iN  -  A  } 

=  e{un_1  Xn_i  I>  <1N  ^  UN  ^ 

=  uN_i  E^N-l  UN  * 


(B9) 

(BIO) 


It  is  easy  to  see  that 

E  ^m-1  ^  =  (WN-1 


h) 


2d 


t 


where 
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